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The master equation for a damped spin well known from the theory of superradiance, is written 
as a finite-difference equation and solved by a WKB-like method. The propagator thus obtained 
looks like the van Vleck propagator of a certain classical Hamiltonian system with one degree of 
freedom. A new interpretation is provided of the temporal broadening of initially sharp probability 
qq"i distributions as the analogue of the spreading of the quantum mechanical wave packet. 
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I. INTRODUCTION 



This is the second of a series of papers concerned with dissipative motion of a large spin which may be realized 
with many identical collectively radiating two-level atoms. The large spin (the "Bloch vector" of quantum optics) has 
conserved length such that its classical state can be described by two angles 6 and if. The classical dynamics is that 
of an overdamped pendulum, 



= sm( 



tan^ = e^tan^, (1.1) 
2 2 ^ ^ 



(P(t) — const . 



o 
in 

• In geodesic jargon we could think of 9 and ip as defining latitude and longitude and speak of creeping motion towards 
I the south pole along a great circle. 
— «. The starting point of the first paper jl] , which we shall refer to as I, was the exact solution of the "superradiance 

^ master equation" (see below) in the form of a Laplace integral |^]. We evaluated that integral in a saddle-point 
r^^*- approximation and derived uniform asymptotics of the dissipative propagator. 
I ' In the present paper we employ a different strategy far less dependent on the specific properties of the problem. We 
O [ observe that in the limit of a large number N of atoms the master equation becomes a finite-difference equation with a 
small step, amenable to solution by an approximation of the WKB type. The propagator solution thus obtained takes 
the form of a van Vleck propagator involving the action of a certain classical Hamiltonian system with one degree 
of freedom. We find the pertinent Hamilton equations to be equivalent to the saddle-point equation of I. The WKB 
approximation entails a new interpretation of the temporal broadening of initially sharp probability distributions as 
the analogue of the spreading of a quantum mechanical wave packet. 
' As in I we shall employ the basis formed by the eigenstates |jm) of and Jz with the respective eigenvalues 
9^ jii + 1) ^^'^ ™- The quantum number j can take on integer or half integer positive values (up to half the number N 
of two-level atoms) and for fixed j the quantum number m runs in unit steps from —j to Denoting the density 
matrix elements by (j, m -\- k\p{T)\i, m — k) ~ Pmi'^) can write the master equation under study as 

= V5m+fe+lffm-fc+lP^+l - (5m - fc^)P™ , (1-2) 

where r is a suitably scaled dimensionless time and g,„ denotes the "rate function" 

5m = j(i + 1) - "'-("^ - 1) • (1-3) 

The particular solution p^(t) satisfying the initial condition p^{t = 0) = 5mn with a certain n is called the 
dissipative propagator and denoted as I?f„„(T). The solution for an arbitrary initial density matrix then is 

As shown in I, the master equation does not couple density matrix elements with different skewness k. In particular, 
the probabilities (k = 0) can be solved for independently of the coherences (fc ^ 0). The elements of the "coherence 
propagator" with fc 7^ are connected by an elementary relation with those of the density propagator [see 



below, (4.1)]. When dealing with the probabilities and their propagator -D^„ we shall drop the superscript k — Q. 
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II. SEMICLASSICAL ASYMPTOTICS 



a. Finite-difference equation — We are interested in the limit of a large number of atoms or, more specifically, of 
a large spin, 



VjU + 1) « j + 1/2 = J » 1 
In that limit it is convenient to introduce a new independent variable /i and its increment A as 

^ = m/J, A = . 



(2.1) 



(2.2) 



In the classical limit /i would become continuous in the range —1 ... 1. For our semiclassical perspective /i remains 
discrete but neighboring values are separated by A. The master equation (l.S) for the densities (case k = 0) becomes 
a finite difference equation. 



dp (/i, t) 
dr 



J 



1 - /i^ - ^A 



p{fi + A,T)- I - / + /lA 



A2 

4 



(2.3) 



b. WKB ansatz — The WKB formalism for finite-difference equations with a small step is well established. The 
general theory has been worked out mostly by Maslov |^ whose notations we shall adhere to. The WKB method for 
ordinary second-order difference equations was extensively used to study the eigenvalues of huge tridiagonal matrices 
occuring in the theory of Rydberg atoms in external field s M . Closer to our topic, the leading (exponential) term in 
the semiclassical solution of master equations of the type ( |2.3| ) was obtained in [|j . We follow the same lines but go a 
step further by also establishing the preexponent, as is indeed necessary to get meaningful results for most quantities 
of interest. 

Let us look for a solution of ( |2.3| ) in a form reminiscent of the WKB wave function in a classically forbidden domain, 

(2.4) 



Here the prefactor A and the "action" S are smooth functions satisfying the initial conditions 

S{p,0) = So{p), A{^i,0)^Aoip) . 



(2.5) 



In our case the absence of the imaginary unit from the exponential does not signal the sojourn of our spin in forbidden 
terrain but simply accounts for the dissipative character of the dynamics in consideration. Incidentally, due to the 
presence of the large parameter J even modest changes of Sq are reflected in wild fluctuations of p(/i, 0); the ansatz 
therefore does not limit our discussion to smooth probability distributions. 

No loss of generality is incurred by assuming the function S (/z, t) independent of J since the prefactor A (fi, t) may 
pick up all dependence on J. We represent the latter by an expansion in powers of A = 



A ifi, t) = {p, t) + (Ai, r) A + A^ + 



(2.6) 



The master equation (2_^) then allows to determine S^A^^^ . . . recursively. We shall need the equations for the action 
and the zero-order prefactor. 



dS_ 
97 



+ (1-A*') 



as 

I — C 3'' 



0, 



d OS 2\ 9 \ , .(0) OS 

efA* (1 - ^'^)— InA'-^'' = c8m 

OT op J 



p. 



(2.7) 
(2.8) 



We shall neglect all higher-order corrections to the zero-order prefactor. 

c. Hamiltonian dynamics — We may consider (2.7) as the Hamilton- Jacobi equation for a classical system with 
one degree of freedom and the Hamiltonian function 



i7(M,p) = (l-M')(l- eP). 



(2.9) 



The canonical equations of motion fi ~ = — (I — /i ) e^, 
resulting "Hamiltonian" trajectories. 



p — —W- = 2pi{l — qP) are easily integrated. In the 
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1 (a + i/)(a-^) 

T = — In — - 

2a [a — ly) [a + ^) 



p = In 



a 



1 -//2 



(2.10) 
(2.11) 



we denote by u the initial coordinate; the name "Hamihonian" is meant to distinguish these solutions from the 
trajectories of the overdamped pendulum (see below). The second integration constant, a, determines the "energy" 
E = H{fi,p) as 



(2.12) 



Rather remarkably, the Hamiltonian trajectory (2.1C) coincides with the saddle-point equation incurred in I when 
examining the asymptotics of the Laplace representation of the propagator. The saddle-point parameter a reappears 
in a new "energetic" role. For later reference we note the nonnegative "speed" 



(2.13) 



A special class of trajectories has zero initial momentum, p(t = 0) = 0, therefore vanishing energy E and a — 1. 
Their Hamiltonian trajectories. 



1 (l + iy)(l- n) 
T = - In ! ) 4-, p(t)^0, 



(2.14) 



are in fact just those of the classical overdamped pendulum (1.1), disguised by fi — cos 6. They involve the canonical 
momentum as conserved with the value zero. 

The semiclassical quantum effects which our Hamiltonian dynamics imparts to the spin through the WKB ansatz 
(2.4) may be seen in the existence of the Hamiltonian trajectories (2.10) not included in the special class ( p. 14 ). 

d. Solution of the Hamilton- Jacobi equation — The familiar relation between canonical momentum and action, 



P 



(2.15) 



implies po (v) = dSo (v) jdv at the initial moment r = 0. Since 5*0 is fixed by the initial density distribution this 
latter equation uniquely associates an initial momentum with the initial coordinate v. One and only one Hamiltonian 
trajectory /i(r; a) thus passes at t = through the initial coordinate v, provided of course that we consider the 
initial probabilities as imposed. Conversely, we can find the initial coordinate v = i^ifJ., t) from which the current 
coordinate fi is reached at time r along the unique Hamiltonian trajectory. 

The action S{fi, t) can now be obtained by integration along the trajectory just discussed, 



S {ij.,t) 



So (t^) + / pdfJ- - Et 



(2.16) 



f—ui^.r) 



We use the explicit form of the Hamiltonian trajectories (2.11) to do the integral. For the sake of later convenience 
we express the resulting action in terms of the auxiliary functions 



a{a, fi, v) = + a) ln(!^ + a) — (/x + a) ln(/i + a) 
— (a — v) ln(a — z/) + (a — /i) ln(a — ji) , 
$(/i, V, t) = [cr(l, /i, v) - a{a, fi, v) + ria^ - 1)] ^^^(^ ,,_^) 



(2.17) 
(2.18) 



(2.19) 



In the definition ( 2. IS ) of the function $ the parameter a must, as indicated in the notation above, be read as a 
function of the initial and final values of the coordinate since these are at present considered as defining a Hamiltonian 
trajectory. We may interprete the function $ as the action accumulated along the Hamiltonian trajectory in question. 
Its derivatives with respect to /i and v accordingly give the final and initial momenta. 
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9$ a? - /i^ 
Oj-l 1 — /i-^ 



dv 1 — 



The Hamiltonian trajectory /i(T; i^, a) can be regarded as the solution with respect to /i of the equation 



dv 



(2.20) 
(2.21) 

(2.22) 



e. WKB prefactor — The expression (2.8) for the prefactor can be simplified using the notion of the full time 



derivative of a function /(/i,T) along the Hamiltonian trajectory /i(r;zy, a), ^'^[^^'^^ — ^'^^^g^'^^'^^ 



- 21 

~ Ot 



since the left hand side in ( p.8[ ) is just the full time derivative dA/dt [see Eg. ( 2. 15 )]. We next introduce the Jacobian 

dfi (r, I') 



Y{t-v, a) 



dv 



y(o,z.) = i, 



and a new exponent B{p,T) to rewrite the prefactor as 

A = 



(2.23) 



(2.24) 



The full time derivative of Y can be transformed to 

dV _ _ d . _ d dH{^i,i 

dr drdu du^ dv dp 



= Y 



dfidp 



dfi^ dp^ 



df, f 






dfidp 


Y exp 


fdS' 





(2.25) 



So equipped we find the simple evolution equation ^ 
the trajectory to give 



B{fl,T) 



-/i for the function B{fi,T) which can be integrated along 
fidT + \nA{v,0) ^ -~^ln "^^ +lnA{v,0). (2.26) 



We thus arrive at the asymptotic solution of the Cauchy problem for our master equation with the initial condition 

(HI), 



1 



piv,0) 



(2.27) 



where v, a are meant as functions of /i and t as explained above. 

/. Narrow versus broad initial distributions — We have in effect constructed the solution (2.27) by the method of 
characteristics. At any rate, the density at a certain point is obtained from the initial density by transport along the 
Hamiltonian trajectory and acquires a factor consisting of an exponential and a prefactor. 

Let us show that the exponential factor is smaller than or at most equal to unity. According to ( 2.20| ) the extremum 
of $ regarded as a function of the final coordinate ^ with v, t fixed occurs when a = 1, i.e. when fi moves like for the 
overdamped pendulum, ^{T;v,a = 1) = ^pend{T]v) as defined by (2.14) or (1.1). At the extremum we have $ = 
and 



d^<^{^i,v,T) 



dfi'^ 

where S is a positive function defined as 



(1-^2) 



(2.28) 



(2.29) 



a— a(/A,^',r) 
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But since $ h as a n egative second derivative its extremum is indeed a m aximum, hence $ < 0. We should mention 
that by using ( 2.2§| ) it is easy to calculate the integral of the density ( ^.27 ) over fi by the saddle-point method and to 
demonstrate that our solution does not violate probability conservation. 

Let us consider two extreme cases of the initial density distribution. First suppose that the initial density po{v) 
is a smooth function, with a gradient of order unity, whence 5o(i^) vanishes. We would have pq — 0, a = 1 which 



means that the characteristic lines are the overdamped-pendulum trajectories ( 2.14 ). The time evolution ( 2.27 ) then 
becomes the fully classical one, 



1 



(2.30) 



The speed of probability transport in this case is obtained by putting a = 1 in ( 2.13 ); since that speed depends on 
the coordinate the initial density distribution will change its form in time, due to its finite spatial size. In particular, 
that change involves broadening or sharpening depending on whether the distribution resides mostly over positive or 
negative values of /i, respectively. 

Now consider the opposite extreme of a narrow initial density, perhaps one almost resembling a delta function. Then 
different parts of the packet will have practically the same coordinate but highly different momenta since po « -j^ . 
The maximum has po = and thus moves along an overdamped-pendulum trajectory. However the parts on the left 
and right slope have, respectively, po > and Pq < 0. Th ey will be transported along the Hamiltonian trajectories 



with a > 1 and a < 1, respectively. But according to ( 2.13 ) that means that the density on the left slope travels in the 
direction of negative fj, faster than the one on the right slope. This will result in a spreading of the initially narrow 
distribution. There is an obvious superficial analogy with the decay of a wavepacket described by the Schrodinger 
equation. Of course, in our dissipative case the exponential factor does not describe any dephasing but rather a 
suppression of probability propagation along trajectories too strongly different from the fully classical overdamped- 
pendulum ones. This puts a brake on the spreading as soon as the packet widens; quantitative estimates of the width 
will be given below in the discussion of the properties of the propagator. 

As an example, in Fig. ^ we demonstrate the time dependence of the probability distribution in the case when 
the system was initially in a pure coherent state of the angular momentum I7) (for properties of such states see 
e.g. [IH^])- The parameter 7 determines the direction of the mean spin vector (7IJI7) as 7 = tan e*"^". We took 
j = 200, 7 = 0.4. Three types of results are presented using: 

1. numerical integration of the master equation ("exact values"); 



2. the WKB solution (|2.27D ; 



3. the fully classical evolution formula ( p.30| ) 



The probability distributions given by the WKB formula coincide with the exact ones with accuracy in the range 
0.4%- 1.7% (accuracy decreases at the later stages of the evolution). The corresponding plots are indistinguishable. 
On the other hand, the fully classical formula correctly places the probability peaks but grossly underestimates their 
broadening with time leading to 20%-150% error in the width and amplitude; this error does not diminish as j 
grows(!) 

g. The semiclassical dissipative propagator — The dissipative propagator establishes a linear relation between 
the initial and final density matrix elements. In the limit of large J the sum in this relation can be replaced by an 
integral; using the classical variables ^, v as arguments it can be written (case fc = 0) 



dvD V, t) p {v, 0) 



(2.31) 



with the function D(/Lt, v, t) related to the matrix Dmn{T) as D [p,, v,t) = J Dmn{T)\j^^j^ n=jv 

To obtain the propagator one has to solve the master equation with t he J -peak as the initial density distribution. 
Strictly s peak ing such an initial condition does not fall into the class (2^) such that our solut ion o f the Cauchy 
problem (2.27) is not directly applicable. It is easy, however, to extract the pro paga tor out of ( 2.27| ) in a slightly 
roundabout way. The semiclassical solution of the dissipative problem in the form (2.4) points to an analogy between 
our master equation for the densities an d a S chrodinger equation in imaginary time. In the spirit of that analogy we 
may consider the function D{p,, v, t) in ( |2.3l[ ) as the van Vleck propagator |^ which must have the structure 



D (/i, V, t) = i?(/x, V, t) 



(2.32) 



with $ (/i, ly, t) the action accumulated along the trajectory. 
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Our task is to establis h th e prefa ctor R. To do so let us substitute the initial density (2.4) and the semiclassical 
propagator in the form ( 2.32 ) into ( 2.31 ) and perform the integration in the saddle-point a pprox imation. The max- 
imum I'* — j/*(/i, r) of the exponent just defines the Hamiltonian trajectories in the form (2.2S). The saddle-point 
integration thus gives 



p(^,t) R{h,v,t) 



27T f [<& {^l,l^,T)+So{u)] 



-1/2 



(2.33) 



where j/*(/x,t) should be substituted for v. Comparing with (2.27) we find the prefactor, 



1/2 



(2.34) 



A simpler form results once we realize from ( ^.22 ) the Jacobian to satisfy 



(2.35) 



dvdfJ. 



For the propagator we thus find 



J 5^$ {fi, v,t) 



271 dvd^ 



1/2 



a? — 



(2.36) 



Compared with the van Vleck propagator for 
an additional square root factor, the origin 
wave functions and density matrices. It is 
that on the classical trajectory (a = 1) this 
V = the inverted classical trajectory 

same factor and combine to the jacobian to 
By expressing the mixed derivative of the 
in ( pgp . 



the Schrodinger equation |g] the preexponential in this expression contains 
of which can be traced to the difference in the no rmali zing condition for 
system specific in as much as a is a solution of (2.11). Note, however, 
factor is just the square root of the classical jacobi an dfj ,~^(^)/dfi, with 
( 2.14 ). One easily verifies that both square roots in ( 2.36 ) give rise to the 
the power one, as it is neccessary to guarantee probability conservation, 
action $ in the preexponcnt in terms of the function S(/i, ly, r) introduced 



a^$(/x, ly, t) 



we arrive at the final form of our semiclassical density propagator matrix Z)„ 



1 



(a2 - ^2)^2^7^ 



(2.37) 



(2.38) 



in which a must be read as the function a = a(fi, v, r) since a Hamiltonian tractory is already uniquely determined 
by specifying the initial and final coordinates. While in general that function can be found only numerically from the 
Hamiltonian trajectories ( 2.1Cl| ), certain limits do yield approximate analytical expressions (see below). Note that a 
is always larger than the larger of 

The behaviour of the propagator solution is largely determined by the action (f>. Two situations deserve special 
mention. First suppose that the initial quantum number n is not close to j. Then we encounter a narrow packet 
centered around the "fully classical" final point UpendiT] i^) where the maximum of $ is located. Close to that 
maximum the propagator can be represented by the Gaussian 



exp 



(M f^pejid^ 
2d2 



(2.39) 



whose width d{h', r) is determined by the second derivative of the action ( ^.28 ). If expressed in terms of the classical 



coordinate fj, the width d — \— J^Lia(Hvend, r)] goes to zero in the classical limit like J such that for many 



purposes the propagator (2.3E) may even be identified with the (5-function. 



6 



A radically different situation is met with when the initial quantum number n is close or equal to j so that the 
initial coordinate v tends to 1 when J — > oo. The action accumulated along any Hamiltonian trajectory is then close 
to zero like 1 — v and so is its second derivative. Indeed, one easily checks what we already saw in I, 



1 - (l - C-2-' j 6, 
l-v ^ {j ~ n)/J 



with 



(2.40) 
(2.41) 



where r' = r — TdassifJ-, i') d enotes a quantum time shift, i.e. the diffe rence between the travel times along the 



Hamiltonian trajectories ( ^.10| ) and the overdamped-pcndulum ones ( ^.14 ). The action then comes out as 



$ « (5 J 1 - 2r' 



,-2r' 



(2.42) 



Obviously, the exponent J$ tends to a finite function of the relative time r' when J tends to infinity as j — n is kept 
fixed. The exponential factor is then only slightly less than unity along a Hamiltonian trajectory. Thus the width of 
the distribution does not tend to zero but rather stays finite as J — > oo. 

h. Comparison of the semiclassical propagators — The semiclassical propagator obtained in I by the saddle-point 
method, 



I^mn ('^) 



1 



(1 - Ai2)V27rJS 



(2.43) 



intriguingly differs from our present (2.3S) in two respects. First, instead of a? — /i^ in the prefactor, ( 2.43| ) contains 



1 — /i^. Second, in (2.43) and in the whole of I we connected the classical coordinate fi to the quantum number m by 



/i = J ^ , which differs from our present definition /i = m/J by the small shift A = J 

Let us show that these two changes in fact cancel each other to leading order in J^^. Consider the form (2.38) of 
our present propagator and substitute /i = /z' + A with /i' the shifted argument (m — 1)/J. If m is not very close 
to n the denominator in the prefactor is not small and one can neglect its change brought about by the replacement 
/Lt — > /i'. As regards the exponent we must exercize greater care: The large factor J obliges us to keep two terms in 
the expansion 



J<^{^l, V, t) = J^{^' + A, i^, t) = J<I>(Ai', V, t) 



In- 



1-^2 



-0(A). 



(2.44) 



The logarithm arising here obviously modifies the prefactor just so as to bri ng a bout t he pr opagator (2.43) of the 
previous paper. Therefore, for J — > oo and m not close to n, the two forms (2.43) and (2.38) of the propagator are 

equivalent. 

If m approaches n the equivalence of the formulae (2.38) and ( 2.43| ) does not hold any more: the accuracy of (IHI) 
is higher. In the extreme cases when n — m or j — n or j ± m are o f order unity or zero |] the uniform approximation 
for the propagator should be used instead of either (2.38) or (2.43) [see I]. 



III. CONCLUSION 



The role of the WKB approximation in quantum mechanics as a bridge to classical mechanics is common knowledge. 
What is not fully appreciated, though, is its potential usefulness for dissipative problems described by master equations. 
The problem of spin damping in superradiance may serve as a goo d ex ample. 

All results of the present paper followed from the WKB ansatz ( ^.4[ ). It led us to a Hamilton- Jacobi equation of a 
classical mechanical system with one degree of freedom. The two canonical variables satisfy two Hamilton equations 
which require two initial or boundary data (like the initial coordinate and the initial momentum or the initial and final 
coordinate) to define a trajectory. This contrasts with the trajectories of the overdamped pendulum model described 
by a single differential equation of first order which are uniquely defined by the initial coordinate only. The extended 
family of Hamiltonian trajectories suggests an intuitive explanation of the quantum broadening of an initially sharp 



^ These cases are similar to the ground state of a quantum system and cannot be adequately described by the WKB method. 
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distributions alien to the traditional classical model: Although all relevant classical trajectrories have the same initial 
coordinate, the initial momenta are different, and so are the final coordinate. 

There is of course a basic difference between the spreading of wave functions and the dissipative propagation of 
probabilities. In the case of unitary evolution the probability amplitudes to arrive from a fixed initial point to 
different final points differ by the phase accumulated along the respective classical trajectory (that is, if wc neglect 
the prefactor). Consequently, propagation along all classical trajectories starting from the same initial point takes 
place with a comparable probability. In our dissipative problem the exponent of the WKB solution is cither zero or 
negative. This means that there is no equality among the Hamiltonian trajectories starting at the same initial point: 
the one which corresponds to vanishing action (the fully classical overdamped-pendulum trajectory) is privileged to 
contribute most. Probability propagation along trajectories significantly removed from that privileged one is effectively 
suppressed. As a result the final width of an initially sharp distribution remains small, unlike the eventually infinite 
dispersion of the wavepacket of a free particle. The exceptional character of the superradiant evolution starting from 
the highest energy levels is also easily understood. In that case the classical action is close to zero on all trajectories. 
Therefore the exponential factor does not limit any more the broadening of the distribution with time. 

The necessity to extend the dynamics of the overdamped-pendulum model was recognized long ago, see the reviews 
1^,0. In a family of trajectories described by two first order differential equations was considered (they were 
introduced as the characteristics in the propagation of quasiprobability distributions) ; the two variables there employed 
did not form a canonical pair, though; one of them was related to our present coordinate /i and the other to the 
transverse components of the spin. 

Our present WKB approach has the merit of being easily extended to other problems in dissipative quantum 
mechanics. Apart from providing an intuitive qualitative picture involving a Hamiltonian equivalent it also provides 
a convenient analytic approximation for the dissipative propagator. In particular we used it to calculate the width of 
the final distribution in the case when the initial direction of the Bloch vector was not close to the z axis. This result 



is of little significance for the theory of superradiance itself. However, the master equation (1.2) becomes increasingly 
important as a model in investigations of quantum chaos in dissipative systems [|l^,|l^ ; disregard of the final width of 
the propagator solutions would lead to erroneous results there. We shall extend the present work to dynamics with 
chaos in a subsequent paper. 

IV. APPENDIX. PROPAGATION OF COHERENCES 

No separate investigation of the coherence propagator (r) is necessary because of the identity proven in I , 

DtSr) = g^n(r) ^^""'''"'^^"'+^'"+" e'^^-/'^ (4.1) 

with 



n _ (i + n)!(j -m)! 
- (j+m)!(j-n)!- ^^'^^ 



It is instructive, however, to consider the changes in our Hamilton-Jacobi formalism necessitated by nonzero k. The 
new quantum number k whose range goes to infinity when j ^ oo must be accompanied with a macroscopic variable 

n^-y (4.3) 



The master equation (L2) written with /i, as arguments now reads 



- J2 (1 - ^2 _ ^2 ^ ^) ^ (5(^2)_ (4 4) 



A Hamilton-Jacobi ansatz 



p{v;fi,T)^A{f^;fi,T) e^^("^^^-) (4.5) 



entails a chain of differential equations for the "action" S and the terms in the expansion of the amplitude A in powers 
of A. We shall examine only the Hamilton-Jacobi equation 
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where F and G denote the auxihary functions 



G{fi)-F{ti)exp(^] = (4.6) 



F(/i) = v/[l-(/i + ^)2][l-(M-^)^], G{f,) = 1 - - rj\ (4.7) 
The previous Hamiltonian becomes extended to 

H{^l,p)^G{^l)~F{^l)eP. (4.8) 
Once more denoting the conserved value of H hy E and introducing the constant a by the relation 



a ^ yJX-E-rf- (4.9) 
we obtain the canonical equation for the coordinate, 

/t = eP = - /i^ , (4.10) 



It coincides with (2.13), and its integration leads to exactly the same trajectories ( 2.10| ) as for the densities. The 



characteristic lines for the coherences propagation are the same as the ones for probability propagation. 

There is one cardinal difference. For A; = we could single out the special trajectories with zero initial momenta. 
Only these were important in the case of smooth initial densities, according to the relation p sa dhi p/ Jdjj.. Since the 
initially vanishing momentum remained zero at r 7^ 0, an initially smooth density remained smooth as long as r was 
not too large (before the system reached its lowest energy level). Smooth density distributions therefore form a closed 
class, and the overdamped-pendulum trajectories are their routes of propagation. 

In the case of coherences it is also possible to select the trajectories corresponding to smooth initial distributions 
or zero inital momenta: the parameter a should be chosen then according to 

= l-if - G{v) + F{u) , (4.11) 

but now these trajectories do not have the physical significance of the fully classical trajectories of the overdamped 
pendulum. The reason simply is that an initially vanishing momentum p will no longer be zero when r 7^ 0. Indeed, 
unless 77 = and a = 1, the momentum 

p = ln ' — = In , ^ 4.12) 

with fi = I^{t) can never be a constant. Remembering the momentum-density connection we conclude that an initially 
smooth distribution of coherences inevitably ceases to be smooth in the course of time. Thus there is no special class 
of characteristics respons ible for transporting smooth distributions of coherences and leaving them smooth, hence no 



elementary relation like ( 2.3C ) 



Our considerations throw light on the important question whether it is admissible to replace the master equation 
by a first-order differential equation using approximations like 

Pin+l ~ Pm + -7: (4-13) 

om 

Such a replacement is justified if the respective elements of the density matrix are and remain smooth functions of 
TO. It follows: 



• In the case of density propagation the above replacement ( 4.13|) is justified provided the initial density distri- 
bution is smooth. 



For coherences the replacement (4.13) is always illegal. 
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FIG. 1. 

Snapshots of the probability distribution p(/i,T) at various times, for the initially pure coherent state (7 = 0.4, j — 
200). The WKB [Eq.(2.27)] and exact results shown by filled contours coincide in the scale of the plot. Dashed 
contours correspond to the classical evolution formula ( ^.30| ) based on the dynamics of the overdamped top. It grossly 
underestimates the width and overestimates the height of the peaks. 
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